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It has recently been pointed out that a substantial amount of e-folds can occur during the waterfall 
regime of hybrid inflation. Moreover, Kodama et.al. have derived analytic approximations for the 
trajectories of the inflaton and of the waterfall fields. Based on these, we derive here the consequences 
for F- and D-term SUSY hybrid inflation: A substantial amount of e-folds may occur in the waterfall 
£SJ ' regime, provided k <C M 2 /Mp, where k is the superpotential coupling, M the scale of symmetry 

breaking and Mp the reduced Planck mass. When this condition is amply fulfilled, a number of 
e-folds much larger than N e ~ 60 can occur in the waterfall regime and the scalar spectral index is 
then given by the expression found by Kodama et.al. n s = 1 — 4/7V c . This value may be increased up 
to unity, if only about N c e-folds occur during the waterfall regime, such that the largest observable 
scale leaves the horizon close to the critical point of hybrid inflation, what can be achieved for 
k ~ 10~ 13 and M w 5 x 10 12 GeV in F-term inflation. Imposing the normalization of the power 
spectrum leads to a lower bound on the scale of symmetry breaking. 

^ '. Among the large variety of inflation models (for a recent review, see e.g. [l|), the hybrid class @, Q is particularly 
promising. In hybrid models, inflation is realized in a false vacuum, along a nearly flat valley of the scalar field 
potential. It ends with a waterfall phase, triggered when the inflaton field reaches a critical value, from which the 
potential in the transverse direction develops a tachyonic instability, forcing the fields to reach one of the global 
minima of the potential. The facts that inflation can be realized at sub-Planckian field values and that it is based on 
■ renormalizable operators only are two attractive features of hybrid models. Moreover, compared to most small field 
models [J], the initial conditions for the fields do not require any extreme fine-tuning because the inflationary valley 
*-0 , is an attractor that can be reached from initial values located outside the valley 

In the usual description of hybrid models @, d, H-Qjl, inflation is assumed to stop nearly instantaneously with 
the onset of the waterfall phase, and the dynamics is restricted to an effective 1-field slow-roll model. Under these 
assumptions, the original version of hybrid inflation exhibits a slightly blue power spectrum and is therefore usually 
considered as ruled out by observations. Moreover, when the Zi symmetry of the potential is broken at the end of 
inflation, domain walls are formed with catastrophic consequences for cosmology. This problem can be solved by 
considering a complex auxiliary field, so that the broken symmetry is t/(l), leading to the formation of cosmic strings. 
In this case, the power spectrum can be in agreement with the CMB data [HI, [l6[ . The fast waterfall phase itself has 
been the object of recent attention, especially to determine the contribution of iso-curvature perturbations [l7l - l22T ] as 
well as the level of non-Gaussianities [23l428j . 

However, it has been pointed out recently that inflation can continue during the waterfall stage for much more than 
60 e-folds [29| , so that the observable perturbation scales exit the Hubble radius during the waterfall. In such a case, 
the power spectrum of scalar perturbations is red, possibly in agreement with CMB observations [Hj-33j, and any 
topological defect formed at the critical point of instability is conveniently stretched outside the observable universe. 

Hybrid inflation can be embedded in various high energy frameworks, like Grand Unified Theories [34|- |36| . string 
cosmology [37r]42j|. extra-dimensions [43|, as well as supersymmetry (SUSY) d, H^-Ell and supergravitv [49l [oOj . 
Supersymmetric models are additionally attractive, because they offer an explanation for the protection of the flatness 
of the inflationary valley against radiative corrections. _F-term [3j,|44| and D-term |45l l49j hybrid models are the most 
well known realization of hybrid inflation in supersymmetry. In these models, radiative corrections lift up the flat 
directions of the potential, giving rise to a red power spectrum of scalar perturbations. In the absence of additional 
non-rcnormalizable corrections, the classic prediction for the scalar spectral index (44| is n s = 0.98. These models have 
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been studied intensively in the effective 1-field slow-roll approach. When the contribution of cosmic strings formed at 
the end of inflation is taken into account, the primordial power spectrum has been shown to be in agreement (even if 
in tension because of large values of the spectral index) with CMB observations [tj [1(3] in some regions of the model 
parameter space. The pred ictions for the spectral index can be lowered when non-renormalizable operators are added 
to the potential [jjl Holrul |48|. l5lj . but the additional parameters limit the predictivity of the models. Moreover, the 
resulting models correspond to hilltop scenarios along the waterfall trajectory, and the initial conditions for the scalar 
fields that lead to phenomenologically viable inflation are less general than for the original models. 

However, the possibility that inflation can continue during the waterfall and affect the observable predictions has 
not yet been explored for the F-term and D-term models. This is the main goal of this paper. We use the method 
of Kodama et al. [30} to integrate the two-field slow-roll dynamics during the waterfall phase and identify that for 
the small coupling regime satisfying k -C M 2 /Mp, where k is the superpotential coupling, M the scale of symmetry 
breaking and Mp the reduced Planck mass, inflation continues for more than 60-e-folds along waterfall trajectories. 

In this case, the standard effective 1-field approach is not valid and the observable predictions are modified. We 
evaluate the amplitude of the power spectrum of adiabatic perturbations, as well as its spectral index. When the 
number of e- folds realized classically during the waterfall is much larger than jV B , the number of e- folds between 
the time when observable modes leave the Hubble radius and the end of inflation, the spectral index is given by 
n s — 1 — 4/iV e . We also calculate the amplitude of the power spectrum, and derive a new constraint on the scale 
of symmetry breaking. In the limit when the observable scales leave the Hubble radius near the critical instability 
point, the spectral index tends to unity, so that it is in principle to find model parameters that accommodate with 
any value of the spectral index in the rage 1 — 4/iV c < n s < 1. 

In Section HH the F- and Z?-term hybrid models are reviewed. A common parametrization of their potential that 
is convenient in order to deal with their dynamics near the critical instability point is introduced. The slow-roll 
parameters are derived and the slow- roll equations of motion are given. In Section [ill! we give a coarse picture of the 
waterfall dynamics that applies when the number of e-folds in the waterfall regime is much larger than 60 and calculate 
the amplitude and the spectral index of the power spectrum of adiabatic perturbations. Section [TV] is dedicated to 
a more precise analysis of the dynamics, and we show, that a spectral index close to unity can arise, provided the 
largest scales, that are observable today, left the horizon during inflation at the beginning of the waterfall phase. In 
SectionfyJ we present estimates for the initial condition for the classical evolution of the waterfall field, that is induced 
by quantum diffusion. Section [VII contains a summary and conclusions. 



II. SUSY HYBRID INFLATION CLOSE TO THE CRITICAL POINT 



A. F-Term Inflation 



The superpotential for F-term inflation is given by 0, |44[ 

W = kS(HH - M 2 ) 



(1) 



where S is a gauge singlet superfield and H (H) are superfields in the (anti-)fundamental representation of SU(A/") 
This gives rise to tree-level scalar potential 

V = k 2 (\HH - M 2 \ 2 + \SH\ 2 + \SH\ 2 ) , 



(2) 



where now S, H and H are complex scalar fields. When S acquires a vacuum expectation value (vev), while (H) = 
(H) = (angle brackets denote the vev), as it is the case along the trajectory that supports hybrid inflation, there 
are Af Dirac fermions of mass kS, N complex scalar fields H + — -^m{H + H) of mass square m 2 + = k 2 (|S'| 2 — M 2 ) 

and AT complex scalar fields H- = -^{H — H) of mass square m 2 ^ — « 2 (|£| 2 + M 2 ). The canonically normalized 

real scalar field a = V2\S\ is the inflaton field, while ip = \/2H + is the waterfall field. Note that the IP-term forces 
\H\ = \H\, implying that the vev of H- is vanishing. 
When (H_) = 0, the tree potential is 



Vb(<7,V) =« 2 M 4 



4M 2 J 



4M 4 



2tf 



The degrees of freedom enumerated above give rise to the one-loop corrections 



Vi 
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log K 



2Q 2 



2Q 2 



2cr 4 log K 



2Q 2 



(3) 



(4) 
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where a c — \f2M is the critical value and Q is an ultraviolet cutoff. SUSY F-term hybrid inflation takes place in the 
potential V = Vq + V%. For a < <r c , the scalar fields are in the waterfall regime, which is the concern of the present 
paper. 

For the dynamics near the critical point a c , the first derivative of the radiatively induced potential, 



dV x {a) 



da 



8tt 2 



; log2, 



(5) 



is of importance. The second derivatives are of order k 4 M 2 / (16tt 2 ) x 0(1). These induce ^-parameters r/ = 
K 2 (Mp/er c ) 2 /(167r 2 ) x 0(1), with Mp the reduced Planck mass. An ^-parameter larger than one violates the slow-roll 
conditions. Therefore, inflation is terminated due to the radiative corrections close to the critical point provided 
k ^> (47T(7 C )/Mp x 0(1). We show below that a substantial amount of e-folds occurs after crossing the critical point 
provided the stronger constraint k <C a 2 /Mp holds. Hence, neglecting the second derivatives is a self-consistent ap- 
proximation in that regime. In the form of the effective potential (U), the 0(1) factor would actually encompass a term 
that logarithmically diverges at the critical point, due to the correction from the massless waterfall field. Moreover, 
when a < <j c , the potential (jlj is ill defined, because of the negative mass-square instability of the waterfall field. This 
is due to the limitations of the method of calculating the effective potential and does not indicate a singular behavior 
in the time-evolution of a. One should expect that the IR-divergence exhibited by the logarithm is regulated by the 
time-evolution of sigma or perhaps the horizon size H . However, even when the logarithm is large, in the regime 
k <C M 2 /Mp, its coefficient is small enough such that we can neglect this effect that would be intricate to deal with 
theoretically. The corrections that lead to the first derivative ([5]) of the effective potential originate from fields with 
positive mass square around the critical point. Therefore Eq. ^ reproduces the slope of the potential at the critical 
point in a reliable manner. 



B. D-Term Inflation 



For this model, the superpotential is [45|, |4£ 



W = kSHH , 



and the D-tcrm is 



D= y -(\H\ 2 -\H\ 2 +m 2 F1 ) 



(6) 
(7) 



The superfields H and H are in the one-dimensional representation of a U(l) gauge group, and mpi is the Fayet- 
Iliopoulos term. The canonically normalized inflaton field is a = y/2\S\ and the waterfall field ip = y/2\H\. For 
(■ijj) = 0, there are two real scalar fields of mass square k 2 ct 2 /2 + g 2 m 2 ?l /A and two of mass square k 2 <j 2 /2 — f/ 2 r7i|j/4. 
When the field a evolves below its critical value 



.Li- 

V2k 



m F i , 



(8) 



the mass square of the waterfall field becomes negative. In addition, there is a Dirac fermion of mass square k 2 <t 2 /2. 
This leads to the tree-level potential 



V = k 2 (\HH 

r 2 



'FI 



2.9 2 



\SH\ 2 - 

^ 2 
2m 2 : 

ff 2 
4k 2 a? 



\SH\ 



l -D 2 



2 

)K 2 a 2 j> 2 
9 «■% 



(9) 



1 i 2 



2k 2 <j, 



When eliminating m-pi in favor of a c , the one-loop potential takes the same form as for the F-term case, Eq. (Q}. 
The ?7-term at the critical point may be estimated as r] = g 2 (Mp / a c ) 2 / (8n 2 ) x 0(1). Therefore we must require that 
g <C 2\f2cr c /Mp, or, equivalently k <C 2mpi/Mp, for inflation not to terminate at the critical point. The comments 
on the IR-divergence of the second derivative and its relevance made above for the .F-term model apply to the present 
case as well. 



Common Parametrization 



Since the F- and D-term models share common features, it is convenient to use the parametrization 



V = A 



1 - a— 



2a 



a 2 ^ 2 



+ Pa 3 a . 



(10) 



The values of the particular parameters can be inferred from the expressions (|3I4I9[) . and they are summarized as well 
in Table |U 
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TABLE I: Parameters to be substituted into the potential (|10[1 in order to obtain the F- and D-term models close to the critical 
point. 



We then follow Ref. [30( in introducing the parametrization 

cr =cr c e^ , 
ip =-0oe x • 



(11a) 
(lib) 



Throughout the slow roll-regime and after the crossing of the critical point, £ < and |£| <C 1, which is consistently 
verified by the explicit solutions. It is useful to note the derivatives 



dV . iaatp 2 . , 
dV 

d 2 v 



da 2 
d 2 V 
dip 2 
d 2 V 



( 4a 2 ip 3 Aaip a 2 - a\ 

[~^T + ~ 2 aT 

4aip 2 



( \2a 2 ^ 2 la a 2 - u\ 
8ao"0 



A 



A 



f4a 2 ip 3 8aip£,\ 



12a 2 ^ 2 8a£ 
7a ^ ~ZT 



dadif) 

and the slow-roll parameters 

1 , _o / 4acr^ z 



£a =-M4 



V 2 



-Mi 



1 , f2 f 4a 2 ip 3 laip a 2 - a 2 



7? CTCT =Mp 



2 - \ (n 

Aaip 2 



-M, 



/4aV 



=Mp 
rj alj , =Af, 



12a^V 4a a 2 - a 2 



2 8aaip 



Aft 



/12aV 8< 



(12a) 
(12b) 
(12c) 
(12d) 
(12e) 

(13a) 
(13b) 
(13c) 
(13d) 
(13e) 



where e x = \M$[(dV/dX)/Vf and t] X y = Af|[9 2 ^/(9A:ar)]/F. The first derivatives enter the slow-roll equations 
of motion, 
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V 



H = TTTT9 ■ (15) 

III. DYNAMICS OF THE WATERFALL: COARSE PICTURE 

In this Section, we determine the field trajectories for the last N c e-folds of inflation that are relevant for CMB 
observations and calculate the scalar power spectrum amplitude and spectral index, in the generic case where inflation 
along the waterfall trajectories lasts for much more than N c e-folds. This regime corresponds to 

aif> 2 /a 2 c « ICI (16) 

and 

4Aa^ 2 > Pot ■ (17) 
It is referred as phase 2(a) in Ref. [3(|. Eqs. (|12a[) and (|12bl) and the slow-roll equations of motion (fH|) yield 



This relation can be integrated, 



dx 2cr2 x 



(18) 



* 2af 2a* [ } 



A sufficient condition for inflation to terminate is the violation of the slow- roll condition |rfy,^| < 1, that occurs for 



T 2 



At this point, rj a ^ = y/2, while 770-0- = (a 2 ) / (8aMp) <C 1. Qualitatively, one may therefore explain the end of inflation 
as a consequence of the classical backreaction via the dimensionless couplings of both evolving fields, a — a c and ip, 
cf. also the discussion in Ref. [20| . Note that the ^-conditions imply that up to this point, no exponential tachyonic 
growth of modes of ip and a has yet occured, such that it is justified to neglect quantum backreaction. 

As a consistency check, we notice that the condition (fT6"|) is met for £ = £ cn d, provided a c /Mp <C 2, such that the 
vev of the inflaton during inflation is sub-Planckian and the effects of (super-)gravity are perturbatively small. At 
this point, we also find r] a ^ — 1 and r\ aiy — a^/(8aMp). This latter term is much smaller than one in the F-term 
model, provided M -C Mp and in the D-term model provided rnpi <C Mp. These conditions coincide with those one 
would impose from demanding (super-)gravity to be in the perturbative regime during inflation. The initial conditions 
chosen for the particular solution (|19[) correspond to a trajectory that trespasses the critical point where £ = and 
if) — 0. At the critical point however, the assumptions (| 161 1 7[) are not valid. Nevertheless, if there are values of £ such 
that £ <C £end and the assumptions (|16I17I) hold on the trajectory (TTO1) . then Eq. (TH)]) corresponds to a trajectory 
evolving from somewhere close to the critical point to £ e nd- 

Using the relation dN = Hdt, where N parametrizes the number of e-folds, and substituting the trajectory p^|) 
into Eq. (|12a[) leads to the equation 



dS, 8M|a 2 



dN 

which can be integrated to 



r, (21) 



a 2 



£ = - (22) 

8aM|(AT end - N + 1) ' v ' 

How many e-folds can occur in this regime? The trajectory p^|) violates the condition (TTT1) . when £ = —^yJN log 2, 
where we set A/" = 1 for the J9-term case. 

Therefore, in order to achieve more than N c e-folds on the trajectory p^|) . the condition 

k « f23l 

2aAf2^A/log2(Af c + 1) V ; 
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must be satisfied. This is one of our main results for supersymmetric hybrid inflation in the waterfall regime. 

In turns out, as described in Section HVl that a substantial amount of e- folds may also occur before the violation of 
condition (fl7|) . In Ref. [3(|, this is referred to as phase 1. The condition for this to happen has the same parametric 
dependence on a, <r c and Mp as for the relation (f2"3")l . but a different coefficient. Before moving to that analysis, 
we derive the predictions for the power spectrum, provided that inflation in the last N c e-folds is described by the 
trajectory (fT9]) . i.e. relation (|23|) is satisfied. 

The instantaneous direction of the evolution of the fields can be parametrized by 

cost? = —=S , (24) 

y<j 2 + ^ 2 

where the dot denotes a derivative with respect to t. We refer to the linear combination of fields in that direction by 
s, which leads to the subscripts that we use in the following. On the trajectory <jTUJ) , cos i3 = l/y/3. The r) parameter 
in the s-direction is 

r] ss = r\aa cos 2 $ + 2r\ a ^, cos ■& sini? + r^y sin 2 $ . (25) 

The parameter r] aa is suppressed by a factor of M 2 /Mp or m FI /Mp, respectively. The other parameters are evaluated 
for ./Vend — N = N e , using Eqs. (|2"Tj) and ([T9|) . For the scalar spectral index, one then finds [3fJ] 

n s = 1 + 2r/ ss = 1 — — . (26) 

iv c 

This result is generic for hybrid inflation in a wide parametric range 30], not only for the SUSY variants. The number 
of e-folds N e is to be evaluated at the fiducial WMAP scale k — 0.002 Mpc -1 . In order to obtain the value for N c , we 
may approximate the energy scale of inflation that is dominated by A. We then obtain [52| 

where ron < A denotes the energy density, below which the Universe is radiation dominated, and where iD ren denotes 
the mean equation of state parameter during the reheating phase (i.e. the end of inflation and the point, after 
which the Universe is radiation-dominated). Its value depends on the details of the tachyonic preheating phase. If 
the tachyonic preheating process is not efficient, u) ro h — and the Universe is matter-dominated due to coherent 
oscillations of inflaton and waterfall field. 

Provided that iso-curvature perturbations do not contribute significantly to the power spectrum, the value for n s 
therefore lies somewhat below its most recently reported central observational value [53L . In order to maximize the 
value of n s , one should therefore assume g rc h = A, corresponding to instantaneous reheating. This can be achieved 
provided the fields H and H have large couplings to the Standard Model sector. Nonetheless, as it stands, the model 
is disfavored at more than 2a even if N ss 60 by CMB observations. However, note that provided less than N c folds 
lie on the trajectory (II9[) , one may expect larger values for n s , because the horizon exit then occurs for vevs at which 
the curvature of the potential is smaller than in the present case. We investigate this in Section IIVI 

Finally, we derive the additional constraint from the amplitude of the power spectrum. The e parameter in s- 
direction is s s = e a + ew, . For the amplitude, we obtain 



A 

V n = 



Jpfcs 

i.e. 



N cnd ~N=N, 



(28) 



k 2 M 2 N i 

P K = £-f for F-term inflation , (29a) 

n 4 M 2 N 4 

Vn = — „ p 9 C for D-term inflation. (29b) 
97r^g^mp I 

In conjunction with the constraint (|23[) . this leads to lower bounds on the symmetry-breaking scales: 

M 2 97V log 2 

-rj2 ^ — TTTH — f° r -F-term inflation , (30a) 

-™P 

m 6 9o 2 lof 2 9 

^ > ^ v& p n for D-term inflation . (30b) 

Alp TT 2 
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These relations together with Eqs. ([29)) constitute another main result for SUSY-hybrid inflation in the waterfall 
regime. 

In FIG. [U we plot the relation between k and M for the F-term model with Af = 1 (notice that the influence 
of the parameter Af is not very significant), that is imposed by the normalization of the amplitude of the power 
spectrum (|29aj) . Moreover, we mark the region in which k is too large (or M is too small) in order to lead to a large 
enough amount of e-folds (|27|) . Lower bounds on the mass parameter M > 10 _6 Mp and on the coupling k > 10 -12 
are deducted. Finally, we have plotted the corresponding energy scale of inflation at the critical point of instability. 
It is found to vary from 10 6 GeV to 10 15 GeV. 



1 




FIG. 1: Relation between k and M for the F-term model with Af = 1 (solid black line), that follows from the amplitude 
of the power spectrum (|29a[l and the number of e-folds The light red region corresponds to the regime of a nearly 

instantaneous waterfall transition, determined with Eq. (123 P . for which no more than a few number of e-folds is realized along 
classical waterfall trajectories. Dashed black lines are the iso-contours of constant energy scale for inflation (V» 1//4 ), respectively 
10 6 GeV, 10 8 GeV, 10 10 GeV, 10 12 GeV and 10 14 GeV from left to right. Blue dotted lines are the iso-contours of constant spectral 
index from Eq. (|2BJ, respectively 0.89, 0.90, 0.91, 0.92, 0.93 from left to right. 



CMB constraints on the three-dimensional parameter space of the D-term model are plotted in FIG. [2] The region 
allowed by the normalization of the power spectrum (I29b[) corresponds to a thin slice of this parameter space. We 
restrict k < 4n, in order to keep the perturbation theory expansion valid. A corresponding two-dimensional diagram 
for various values of g is provided in FIG[3J For 0.1 Mp < mFi < Mp, the spectral index can be in agreement (but in 
strong tension) with the WMAP constraints, provided a coupling to fermions of the order of unity and n ~ 1CU 4 . 

Let us notice also that the energy scale of inflation for the D-term model can be as low as a few TeV and the 
number of e-folds during the waterfall about N c , provided g « 10~ 20 , n « 10~ 16 and mFi ~ 10 _8 Mp. This extreme 
case is of particular interest since the D-term model can then provide a mechanism for the recent acceleration of the 
Universe expansion (55j . 



IV. DYNAMICS OF THE WATERFALL: REFINED PICTURE 



We apply now some more of the details that are derived in Ref. [30( to SUSY hybrid inflation. In particular, we 
consider initial conditions that satisfy relation (TT5]) but that turn around the bound (fT7| . such that 

4Acn/> 2 < Pa & c . (31) 

In Ref. [3(|, this is referred to as phase 1(a). Use of the slow-roll equations of motions (ITil) and Eqs. (|12a|) and (|L2b|) 
leads to the differential relation 

* - M. (32) 
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FIG. 2: Constraints on the 3D parameter space of D-term model from the normalization of the power spectrum (|29b[l (blue 
area). The red region corresponds to the regime of nearly instantaneous waterfall, determined with Eq. (I23|l . 




1(T 6 1(T 5 1(T 4 0.001 0.01 0.1 1 

m Fl /M P 

FIG. 3: Relation between k and mpi for various values of g (from top to bottom: 1, 10 -3 , 10 -6 , 10 -9 , 10 -12 , 10~ 15 ) in the 
D-term model from the normalization of the power spectrum (|29bl) (blue lines). Above the red lines, that are determined by 
Eq. (|23[) . the waterfall transition proceeds nearly instantaneously. Black dashed lines and blue dotted lines are iso-contours of 
constant energy density and spectral index values, as in FIG[T] 



what can be integrated to 



(33) 
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Here, £1 and xi should be determined by appropriate boundary conditions. The point on this trajectory, where the 
condition (f3"Tj) is violated, is denoted by £2 and Xi- When £ 2 ^ £1 and X2 3> Xli on e may approximate £1 w 0, Xi 1=3 0j 
which is what we assume in the following. Using these approximations, we obtain 



1, _K 

"2 g 4aA^ 



=0 log 7^2 > (34a) 



£=5S»- (34t) 
The number of e-folds before reaching (£21X2) is 

" a =2^^' (35a) 
M 2 /~2y7~ 

N 2 = 2ttJ-^- for F-term inflation , (35b) 

kM 2 V log 2 v ; 

^2 =^|-27r. /-^- for D-term inflation . (35c) 

KAf 2 yios 2 

Therefore, also in the regime where relation (f3"Tj) is valid, a substantial amount of c-folds may occur in the waterfall 
regime provided K <C M 2 /Mp or k <C m|j/Mp, respectively. 

Substituting the trajectory (|33p into Eqs. (|12all4[) and making use of the relation (I5T1) . one finds 

M 2 

e(iV) = -iV^/3 ( 7 c 2 ) (36) 

such that the fields are close to the critical point for N = 0. 

We check whether before reaching £2, the condition (fT5|) may be violated. The inequality (TT5|) holds for all points 
on the trajectory (f3"3")l before reaching (£27X2)7 provided that 

X2 »^"^ , ( 37a ) 
K 2 

y 2 > ^log2 for F-term inflation , (37b) 

A 128tt 2 6 ' v ; 

<? 4 

X2 o log 2 for £>-term inflation . (37c) 

2oott z k z 

Since xi is given the by logarithm in Eq. (|34aj) and k -C 1, it is immediately clear that above condition holds for the 
F-term case (barring the choice of large values for V>o)- For D-term inflation, Eq. (|35c|) and relation (|37cl) combine to 

- » ~^hr , (38) 

which is satisfied because we already observe the stronger constraint a c -C Mp n/g mpi/(\/2Mp). In conclusion, 
the condition (fT6)) is fulfilled for both, F- and Z?-term inflation at all times during phase 1(a), and when combining 
this with the results of Section HTH it is fulfilled at all times during the slow-roll regime as well. 

After the transition from phase 1 to phase 2 [the point (£21X2)]) condition (|3"Tj) no longer holds and is replaced 
by (|17l) . The fields £ and x satisfy the differential relation (|18[) . 

Compared to the trajectory (|19l) . a solution can be determined that takes accurate account of the boundary 
conditions that arise at the end of phase 1 [3(| : 



-{d ~c)f{N)+c' +c cl 
{cJ - c)f(N) +c' + c 

f(N) = e lec ' M 'pVW(N-N 2 ) ^ (39b) 

c= v ^ 2 72, (39c) 



r= v <?-~, (39d) 
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At late times, the solution (|39al) approaches the approximate form (|2Tj) . Notice that the initial conditions (&,X2) 
depend on the initial condition for phase 1 through while the late-time behavior (|2ip is independent of these. 
We can substitute £ cn d, Eq. (J20J) , in the solution (|39a|) . invert it and obtain the number of e-folds in phase 2 [30l |: 



A^nd -N 2 = 



1 



a/3 VCond+?2 C - C ' 



(40) 



which corresponds to an improved version of the estimate (|23[) . 

0.0000 

-0.0002 - 

-0.0004 - 

-0.0006 



-0.0008 



-0.0010 




FIG. 4: Field trajectory in the space for F-term inflation, with k — 0.001 and M = O.lMp. The plain blue curve is 

obtained from Eqs. ()22I19|) It is in agreement with the red dashed curve, obtained by integrating numerically the exact classical 
dynamics. 



It would be straightforward now to derive a general expression for the power spectrum in a situation where -/Vend — 
N2 < N c , such that the horizon exit of the largest observable scales occurs in phase 1(a). The result is however 
somewhat complicated and shows no advantage over a numerical evaluation, while having the disadvantage of incurring 
an inaccuracy due to the matching of phases 1(a) and 2(a). There is however an interesting limiting case, where simple 
estimates can be obtained: While in Section IIII1 we have discussed the situation where the number of e-folds during 
the waterfall regime is much larger than N c , i.e. N cn a ^> iV e , we now consider the situation where just enough e-folds 
occur in the waterfall regime, i.e. N cn a > N e and N en ^ — A" c <C N c . When e„ ^> e^,, as it is the case sufficiently close 
to the critical point, the amplitude of the power spectrum in phase 1(a) is 



A 



A 3 



2A-k 2 M 2 e 



(41) 



The prediction for the scalar spectral index now depends on the initial condition for ijj, i.e. the choice of ipo in the 
present parametrization. However, as explained in Sec.[V]below, the auxiliary field dynamics at the critical instability 
point is not classical but dominated by quantum diffusion effects. The classical regime is nevertheless reached quickly 
and our initial value i/jq must be seen as the value that would take the auxiliary field if the classical trajectory was 
traced back up to the critical instability point. 
Note moreover that when a ip, 



a .,4 32a 2 V 2 



costfwl-iV^M, 



(42) 



When we impose that the exit scale is close to the critical point, the superpotential coupling must be of order 
k ~ M 2 /Mp orK~ m|j/Mp. Close to the critical point, when assuming that ip is given by a value close to its lower 
bound (f53"| . it then turns out that the effective ^-parameter (f2"5| is suppressed by M 2 /Mp or m Fl /Mp. Therefore, 
the scalar spectral index approaches values very close to unity when the horizon exit of the largest observable scale 
occurs close to the critical point of SUSY hybrid inflation. It should therefore be possible to find parameters, that 
produce all possible values for the scalar spectral index n s between the values (|26l) and one. The parameters for which 
the horizon exit occurs close to the critical point can be estimated from FIGS. [TJ [2] and [3] from the intersection of the 
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exclusion region, for which less than iV c e-folds occur below the critical point and the relation between re and M or 
inpi, respectively. In particular, for F-term inflation, we can estimate that re s» 10~ 13 and M » 10 12 GeV, in order 
to yield a value of n s close to its presently observed central value. A precise determination of this point in parameter 
space and a quantitative analysis of how much tuning is required, such that the spectral index falls within the allowed 
range, will be subject of a future study. 

Calculating the possible values of re, M or mpi for a given n s using the present methods does not appear to be 
possible in a simple analytic way. For example, Eq. f|40[) cannot be solved for re in terms of elementary functions. Due 
to the simple nature of the slow-roll equations, a numerical study should however be feasible. 



V. FROM QUANTUM DIFFUSION TO CLASSICAL DYNAMICS 

The statistical distribution of the initial auxiliary field values ipo can be evaluated by studying the quantum diffusion 
near the critical instability point [2i| The coarse-grained auxiliary field can be described by a Klein-Gordon 
equation to which a random noise field X(t) is added (5?j]. This term acts as a classical stochastic source term. In the 
slow-roll approximation, the evolution is given by the first order Langevin equation 

1 dF H 3 / 2 

^ + 3i7d^~ A W' ™ 



which can be rewritten by using Eq. (|12b[) 



2 / a 2 



■ H 3 / 2 , , x 4cW>M v , 

2vr a 2 V a 2 ) 

The two-point correlation function of the noise field obeys 

<A(f)>=0, (\(t)\(t'))=S(t-t'). (45) 

In the limit of H constant (this approximation is valid at the critical instability point when the expansion is governed 
by the evolution of a in the false vacuum), this equation can be integrated exactly. Under a convenient change of 
variable 

x =^ =exp (-2N¥zfi<A , (46) 



A 



one gets 



dip _ H 1 ' 2 w ^ 4aipM$(l - x) 



dx inrx^ ^ 2a 2 rx ' ^ ^ 

with r = Mpf3a 2 /A This equation has an exact solution 
tp(x) — C exp (C2X — C2 lax) 

(48) 



C\ exp (C2X — C2 In x) x J exp [—C2X + C2 In x') X(x')dx' 



where C\ = H 1 / 2 /(Anr), C2 = 2a/(a 2 r) and C is a constant of integration. The variance of the auxiliary field 
distribution is then obtained by taking the two point correlation function of ip(x). By assuming an initial delta 
distribution for ip at cr 3> <r c , one obtains 



H 



2 r 



exp(x) 



T(a,ax) , (49) 



where we have defined a = 4aMp/(a 2 r) = 4aA/(/3er^) and where T is the upper incomplete gamma function. Near 
the instability, x ~ 1 and one thus has 

u"2 

(^ 2 (xcl))^^-(^) a r(a,a). (50) 
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By using recurrence relations as well as the asymptotic behavior of the T function, one can hnd 



(e \ u Pk ] 

— ) xr(«,«)~(/- — = when u — > oo , (51) 
uJ y 2 \Ju 

such that 

(0 2 (x~l)>~ ^~^ = f'ft (52) 

At the critical instability point, the average value of -0 over the whole Universe is zero, and Eq. (|52"|) describes the 
statistical distribution of the held around zero. However, over a small patch that will contain our observable Universe, 
the average value is non zero and increases statistically with time due to the second term of Eq. (|44|) . But the variance 
(ip 2 ) in this patch is still given by Eq. (|52|) . After some e- folds of inflation, the classical regime is reached and the 
classical evolution of rp proceeds faster than quantum diffusion, H^ip > H . By using Eqs. (|12bll4l31l36[) . one finds 
that this happens when 

A 3/2 

1> > —r= ? • (53) 

We should therefore see the classical dynamics of tp during the waterfall as emerging in a patch of the Universe where 
the quantum diffusion was previously dominating. But the time when inflation takes place, the dynamics of the held 
cr remains classical. In particular, when combining the estimates (|52l) and (|53")) . one may see that typically, N <C 1 
when the classical evolution begins: for D-term inflation, N ~ nM 6 /Mp 6 and for D-term inflation, N ~ gm Fl / Mp 6 . 

Besides the problem of the quantum diffusion of the auxiliary field, one must also take care that the inflaton itself is 
classical. For the original hybrid model, the regime dominated by the quantum stochastic fluctuations of the inflaton 
has been studied in Ref. [56| and leads to a strong reduction of the number of e-folds realized during the waterfall. 
This argument imposes the additional condition |dcr/diV| ^> H/(2ir). During the phase 1(a), the classical evolution 
of a is governed by (1361) , so that this condition can be rewritten 

» 1 . (54) 

A 3 

It is satisfied provided KMp/M 3 ^> 1 for the D-term model, and nMp/m Fl ^> 1 for the D-term model. It therefore 
appears for the D-term model that in the range O.lAfp < mpi < Mp, g ~ 0(1) and k 10~ 4 that was found to be in 
agreement with CMB observations, the quantum effects of a during the phase la can be important. However, in that 
particular case, since the inflaton is driven by the second term of (|12aj) during the last 60 e-folds of inflation (phase 
2a), we argue that its dynamics is classical during this phase such that observable predictions are not affected by the 
quantum stochastic effects of a at the critical instability point. 



VI. SUMMARY AND CONCLUSIONS 



In the present work, we have explored the prospects of SUSY F- and D-term models of accounting for the observed 
normalization and spectral index of the primordial perturbation power spectrum. We have focused on the parametric 
regime where all scales that are observable today have left the horizon during the waterfall stage. It is particularly 
interesting to confront these scenarios with observational data, because they only rely on the scale of symmetry 
breaking and the superpotential coupling for D-term inflation and, in addition, the gauge coupling for D-term inflation. 
When the scale of symmetry breaking is small compared to the Planck scale, the influence of non-renormalizable 
operators, that is expected within the supergravity completion of these models, is suppressed. Besides, the trajectories 
of the scalar fields are attracted to the valley of the hybrid inflation potential, what makes these models rather 
predictive. 

In order to derive our results, we use the analytical methods that have been introduced in Ref. [30| . A more accurate 
numerical study will be subject of future work. The main conclusions, that we presently achieve, are as follows: 

• Inflation proceeds in the waterfall regime, provided k <C M 2 /Mp 2 or k <C m Fl /Mp 2 , respectively. More accurate 
relations are given by Eqs. (|23l) and (HU1) . 

• When the number of e-folds of inflation in the waterfall regime is much larger than 60, the dynamics can be well 
approximated as in Section [TTT1 There is no restriction on the scale of symmetry breaking, but the normalization 
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of the power spectrum imposes a relation with the superpotential coupling k, cf. Eqs. (|29ll4Tj) and FIGs. [Tl [2l [3] 
The spectral index n s , Eq. (f2l)f (cf. also Ref. [30j] ) , then takes values below its present central observational value. 
It can be in agreement (but in strong tension) with WMAP only for the D-term model with O.IMp < ttifi < Mp, 
g ~ 0(1) and k « 10~ 4 . 

• It is therefore interesting to study how larger values of n s can be achieved. In Section IIV1 we have shown that 
provided the largest observed scales leave the horizon close to the critical point, the deviation of n s from unity 
is suppressed as M 2 /Mp 2 or m FI /Mp , respectively. The point, where just enough e-folds of inflation in the 
waterfall regime occur and the power spectrum is normalized in accordance with the observed values can be 
inferred from FIGs. Q] O |3l In particular, for F-term inflation, we can estimate that this situation occurs for 
K w 1CT 12 and M « 5 x 10 12 GeV. 

When the observed limits on n s further tighten around its presently observed value, it is in order to further study 
the latter possibility. For the purpose of determining the parameters more accurately than by the order-of-magnitude 
estimate in the present work, it will be necessary to perform a numerical study instead of the present analytic 
approximations. In particular, this is necessary because of the inaccuracies in the matching between the phases 1 
and 2, as explained in Section ITVl and because we did not consider the possible contribution of iso-curvature modes. 
An interesting question will be how the observational uncertainty in n s will translate into allowed ranges of k and 
M or ttifi, because this will quantify the amount of parametric tuning that is required for SUSY hybrid inflation to 
match observations. 

Finally, we have mentioned for the D-term model an extreme case in which the energy scale of inflation is only a 
few TeV, so that the model can provide a mechanism to the current cosmic acceleration of the expansion, accordingly 
to Ref. [55|. This occurs for g w 1CT 20 , k « KT 16 and mm ~ KT 8 Af P . 



We conclude that SUSY hybrid models will remain interesting proposals in order to explain the observations of 
the primordial power spectrum, even when it is further confirmed that n s is substantially below 0.98. Even without 
effects from additional non-renormalizable operators, smaller values of n s are predicted provided the largest observed 
scales left the horizon during the waterfall regime. Note that also in this situation, the appealing features of the 
model, which are the dependence on a small number (2 or 3) of renormalizable operators only, the natural emergence 
of inflationary field configurations due to the attractor property of the potential and finally, the motivation from 
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